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First Passage and Cooperativity of Queuing Kinetics 
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We model the kinetics of ligand-receptor systems, where multiple ligands may bind and unbind to 
the receptor, either randomly or in a specific order. Equilibrium occupation and first occurrence 
of complete filling of the receptor are determined and compared. At equilibrium, receptors that 
bind ligands sequentially are more likely to be saturated than those that bind in random order. 
Surprisingly however, for low cooperativity, the random process first reaches full occupancy faster 
than the sequential one. This is true except near a critical binding energy where a 'kinetic trap' 
arises and the random process dramatically slows down when the number of binding sites N > 8. 
These results demonstrate the subtle interplay between cooperativity and sequentiality for a wide 
class of kinetic phenomena, including chemical binding, nucleation, and assembly line strategies. 

PACS numbers: 02.50.Ey, 05.20.Dd, 05.60.-k 
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Cooperativity plays a key role in determining the equi- 
librium properties of queuing systems such as ligand- 
receptor binding, nucleation, melting of a-helices, cou- 
pled chemical reactions, and assembly lines. Biophysi- 
cal examples include O2 or CO binding to hemoglobin 
and myoglobin P, 0,0] and binding during cell signaling 
and morphogenesis 0, Q . For a local bulk ligand con- 
centration, the associated receptors will typically have a 
fraction of sites filled. The kinetics of queuing in these 
processes can exhibit diverse and rich behavior. A recep- 
tor may need to have a critical number of bound ligands 
before it can signal the next biochemical step. Thus, it is 
important to know not only the equilibrium ligand occu- 
pancy, but also the mean time to first reach this critical 
occupancy, as a function of local ligand concentration 
and binding strength. Similarly, in nucleation processes 
such as cv-helix formation or melting, local helix turns can 
form randomly or sequentially. The first time a complete 
helix forms (or melts) will be an important ingredient 
in protein folding models 0. First passage times also 
define extinction and fixation in birth-death processes 
0, @ • Equilibrium distributions and first passage times 
also arise in applications of queuing, where, for example, 
average computer loads and the first time that demand 
exceeds capacity should be distinguished HJ. 

In this Letter, we formulate and use a kinetic chain 
model to highlight subtleties of ligand adsorption and 
desorption, queuing, and cooperativity Our model is 
presented in the language of ligand binding to a single 
receptor with N active sites of which < n < N are 
occupied by ligands at any given time. The order of the 
binding can be imposed in two limiting ways. As shown 
in Fig.^ the addition of each successive ligand can in- 
fluence one other specific site and allow the next ligand 
to bind to, or unbind from, that site only (a case we will 
denote by the index a = 0). Alternatively, the allosteric 
effect (from e.g. a large scale conformational change) can 
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FIG. 1: Two limiting models of multiple ligand binding. In 
both cases, N — 6, and n = 3. (a) Sequential ligand bind- 
ing/unbinding. At any given time, only one specific site can 
accept ligand adsorption and only one ligand can desorb. 
Upon defining the kinetics, this case will be denoted a — 0. 
(b) Random ligand binding/unbinding. An additional ligand 
can bind to any one of the open sites. Similarly, any one of 
the bound ligands can spontaneously desorb. This case will 
be denoted a = 1. 



be spread equally to all remaining sites. The next ligand 
can bind to any one of these remaining open sites (a case 
we will denote by the index a = 1). Here, all bound 
ligand molecules are equally likely to spontaneously des- 
orb. We do not consider mixed processes in which the 
binding order is sequential and the unbinding random, 
or vice versa. Thus, in our model, ligand binding occurs 
in a totally sequential manner as in Fig.^i, or randomly 
as shown in Fig.f^D. We show how the binding order 
plays a crucial role in determining both equilibrium and 
kinetic properties: sequential ordering generally implies 
higher occupancy, while random ordering is associated 
with a shorter mean first passage to saturation. We also 
discover an intriguing regime of kinetic slowdown where 
different queuing rules result in dramatically different be- 
haviors. 

We define P(n, t) as the probability that the receptor 
has n bound ligands at time t, given that it had hq at 



time to. The evolution of the occupancy state n can be 
mapped onto a one-dimensional random walk with a mas- 
ter equation given by d t P(n, t) = q n+ iP(n + 1, t) — (q n + 
k n )P(n,t) + fc,,_iF(n — l,t). Here, k n and q n represent 
the ligand adsorption and desorption rates, respectively, 
when there are already n ligands bound to the receptor. 
The probability density current from state n + 1 to state 
n is J(n,t) = q n+ iP(n + l,t) — k n P(n,t). The master 
equation, d t P{n, t) = J(n, t) — J(n — 1, i), is solved under 
the constraint Y^iLo-Pi = 1- 

Different functional forms of k n and q„ distinguish the 
sequential model from the random one. The number of 
ways for a ligand to bind an n state receptor is propor- 
tional to the number of accessible binding sites. For se- 
quential kinetics, at any given time, only one site is ac- 
cessible for binding the next ligand, whereas for random 
binding, all N — n empty sites are accessible. The num- 
ber of accessible sites can thus be succinctly expressed 
as (N — n) a , where the index a = corresponds to the 
sequential model, and a = 1 corresponds to the random 
one. Similarly, the number ways to detach a particle 
from an occupation state n + 1 is (n + l) a . If conforma- 
tional changes of the receptor molecule reach local ther- 
modynamic equilibrium within the time scales required 
for ligand binding and unbinding, the rates obey detailed 
balance: 
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FIG. 2: Equilibrium occupation probabilities P eq (n) for N = 
10 and 7 = as a function of ligand loading n and various 
affinities x. (a) Sequential binding, (b) Random binding. 
Only the random model exhibits a probability maximum at 
intermediate n. This results from enhanced opposing drifts 
at small and large occupancies. The distributions for 7 > 1 
are qualitatively similar to those in (a), (c) The equilibrium 
probability of full occupancy (n = N — 10) for both models 
and 7 = as a function of x. Full occupancy is more likely 
for the sequential case than for the random one. (d) The 
sequential and random case MFPT at TV = 10 and 7 = 0. 
Note the logarithmic scale. 



The proportionality constant z w v[L] is dimcnsionlcss 
and depends on the bulk ligand concentration [L] and on 
the capture volume per binding site v. Typical values of 
[L] = ljiM and v = Inm 3 give z ~ 10 -6 . The quantity 
— AG n = G n — G n +i is the free energy change, in units 
of thermal energy fe^T, upon detachment of one ligand 
from state n + 1 to state n. Equivalently, AG n can be 
interpreted as the free energy change due to ligand ad- 
dition from state n to state n + 1. Cooperativity can 
be defined as the additional proclivity for the (n + l) st 
ligand to bind as n increases. A simple model for the 
free energy of binding is AG n = —Eo(n + l) 7 where Eq is 
the nonncgative ligand binding energy in units of ksT. 
The nonncgative parameter 7 controls the cooperativity 
of successive binding. For 7 = 0, AG n = — Eo is in- 
dependent of n and the binding is noncooperative. If 
7 > 0, AG n decreases with n, hence a positive cooper- 
ative effect. Large values of 7 represent strong cooper- 
ativity. The equilibrium probability distribution P eq (n) 
is derived by imposing J(n, t) = and by using reflect- 
ing boundary conditions for the empty and full states, 
fc_i,<7o = and gjv+i, few = 0, respectively ^(|- Flpon 
using Eq.^ 
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For the noncooperative case, the occupation probabil- 
ities can be expressed in terms of a binding affinity 
x = ze £ °. In Figs.Ek-b wc show the noncooperative 
P eq (n) for N = 10 and various x. Only the random case 
allows an intermediate maximum of P eq (n) to develop as 
x is increased. From Eq.Ql it can be seen that for small 
n values the ratio k n /q n+ i is larger for the random case 
than for the sequential one. Conversely, the same ratio 
is smaller for large n. The random process thus induces 
an effective drift towards intermediate n values. This is 
schematically shown in Fig.|3| We can rewrite Eq.^ in 
the form k n /q n+ i = cxp(— AG n ), where 
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FIG. 3: Kinetic steps for x ~ 1 and 7 = 0. (a) Sequen- 
tial process. All successive steps incrementing the number of 
bound receptors have the same rate k, while all decrement 
steps have the same rates q. (b) Random process. There 
are more ways to attach a ligand when n is small. Similarly, 
there are more ways of removing receptors when there are 
many bound receptors. The effective rates push the system 
from both directions toward an intermediate occupation level. 
Given an initial condition at n = 0, the first passage statis- 
tics to n = N will be controlled more by the effective forward 
rates k n than the effective backward rates q n . Therefore, the 
larger k n tending to load the receptor outweigh the effects of 
the larger q n tending for desorption. The overall MFPT is 
thus larger for in the random case than in the sequential one. 



is the effective free energy change that includes the en- 
tropy associated with ligand binding and unbinding. Un- 
like G„, the total effective free energy as a function of n, 
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can have a minimum with respect to n if a — 1. Thus, 
a maximum in P eq (n) can occur only in the random 
case. Equilibrium distributions and occupancy probabil- 
ities have been widely used in approximations of kinetics 
of biochemical reactions |icj . One commonly used metric 
is the filling fraction f eq (x;a,j) = N^ 1 J2n=i n Peq(n). 
We can compute f eq under either the sequential or ran- 
dom assumptions: 
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respectively. These two expressions coincide only in the 
infinite affinity limit x — > oo, signifying that the process 
no longer depends on binding order. Expressions such as 
f eq have been used in kinetic equations describing bind- 
| to hemoglobin/myoglobin EIQ, regulatory networks 
viral infection dynamics |l2^ , and cell signaling |l3l | . 
However, using equilibrium expressions in kinetic equa- 
tions assumes separation of time scales or near steady- 
states. Moreover, filling fractions do not provide infor- 
mation about typical loading time scales. The queuing 



noncquilibrium properties arc more telling, in this con- 
text, than their equilibrium counterparts and surprising 
effects may arise. In particular, faster complete occu- 
pation of a receptor can occur even though the mean 
occupancy is lower at equilibrium. The mean first pas- 
sage time (MFPT) to a target filling can be evaluated 
by using standard techniques Ql- F° r a receptor in the 
unoccupied n = state at t = the MFPT to reach an 
n* filling is flH | 
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To evaluate Eq. Q we need independent expressions for 
k n and q n . We let k n = vv[L](N — n) a and q n+1 = 
v (n + l) a exp [AG„], where v is the frequency of reac- 
tion attempts. We choose k n to be independent of AG n 
because for many systems, such as antibody-antigen com- 
plexes , cooperative effects are felt only in the detach- 
ment process. In these cases, the strength of a ligand- 
receptor chemical bond depends on the total number of 
bound ligands present. These existing ligands typically 
do not affect the kinetics of attachment. For the nonco- 
operative case (7 = 0), the MFPT to full occupancy in 
units of v~ x is jl7j ]: 
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Logarithms of the MFPT from n — to complete filling 
n* = N are plotted in Fig.|2jl for 7 = 0. A compari- 
son with Fig.[2i shows an unexpected effect. Although 
the noncoopcrative random process yields a lower equi- 
librium ligand population than the sequential process, its 
MFPT is shorter. This is due to an overall effective in- 
crease in the local diffusivity as evident from the form 
k n oc (N — n) a . However, near the half-filling affinity 
x = 1, corresponding to f eq = 1/2, the MFPT of the 
random case is greater than that of the sequential one. 
This is emphasized in Fig.^Ji, where we plot the ratio 
of random to sequential MFPT as a function of binding 
energy. A striking peak develops for x = 1 provided the 
receptor has at least A^ c = 8 binding sites. Indeed, as 
N > N c increases, the random MFPT increases dramat- 
ically. Why does this happen? From Fig-El we see that 
for the random case near x = 1, the backward rates q n +i 
balance the forward rates k n . This results in a vanish- 
ing net drift at intermediate occupations. The effective 
potential G n develops a minimum that for x = 1 is ex- 
actly located at N/2. The opposing drifts bias the filling 
to intermediate values forming a kinetic trap and delay- 
ing first saturation. The MFPT increases significantly 
due to the time required to escape this trap, despite the 
increase in effective diffusivity near the potential mini- 
mum. When N > 8, the trapping is strong enough that 
the MFPT of the random case is greater than that of the 



sequential process. This effect is disrupted as x deviates 
from unity as the minimum in G n disappears. Note that 
this behavior occurs only in the absence of cooperative 
effects: for e. g., 7 = 1 the equilibrium occupancies fol- 
low the same trends as in Fig.|2t and the random system 
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FIG. 4: Mean first passage time ratios T(N; 1)/T(JV; 0). (a) 
For 7 = 0, random binding (a — 1) is faster except for a 
window of binding energies near x = ze e ° — 1. For clarity 
we plot T(N;1)/T(N;0) as a function of e for z = 10" 6 . 
The critical value x = 1 occurs at £0 = 13.8. The peak 
scales as ~ 2 N /N 3 . (b) For 7 = 1 (and in fact 7 > 1), the 
random process always reaches full occupancy faster than the 
sequential one. 



To explicitly see a size-dependent trap, consider the 
MFPT, T trap (M; a), out of a band of sites of length 2M, 
centered about n = N/2. For x = 1, 7 = 0, and an 
initial condition of n = N/2, the logarithm of the ratio 
of the random T trap (M; 1) to the sequential T trap (M; 0) is 
plotted in FigJSJ The trapping effect is evident from the 
fact that for large enough M the random T trap (M; 1) can 
easily be e 10 times larger than the sequential T tra p(M; 0). 

Our analysis reveals the conspicuous effects that co- 
operativity and sequence have on binding kinetics, es- 
pecially on the MFPT to saturation. The key result is 
that equilibrium and dynamic measurements offer differ- 
ent answers as to whether the sequential or the random 
order is more efficient in ligand saturation. At equilib- 
rium, sequential processes are more likely than random 
ones to saturate the ligand. For dynamic properties, 
such as MFPTs, random processes are faster. This is 
true except for large noncoopcrativc systems (N > 8) 
with a binding affinity near x = 1, where random pro- 
cesses dramatically slow down. One can also consider 
the reverse process of emptying a completely filled re- 
ceptor. The clearance time can be obtained by using 
particle- hole symmetry and the replacements k*/q* +1 — 
[(N — p)/(p+ l)] a exp(AGp) / z where p is the number of 
holes present. 

The rate parameters used in our study can also be 
reinterpreted as the inverse of the mean servicing time 
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FIG. 5: Logarithm of the ratio T t ra P {M; l)/T trap (M; 0). 
Here, x = 1, 7 = and the initial occupancy is at n = N/2. 
Negative values of the ratio imply that escaping a 2M kinetic 
trap is faster for the random process than for the sequential 
one. This occurs only for small enough traps centered on N/2 
where the random dynamics is faster. For large enough val- 
ues of trap length 2M the effects described in Fig.^take place 
and the exit time of the random case, Ti rap (M; 1), is much 
larger than the sequential one, T tra p(M;0). 



in a customer service queue. Our results indicate a rich 
set of outcomes depending on queue size, cooperativity, 
and sequentiality. Levels of "service discipline," or how 
to order customer service @ , can be implemented in dif- 
ferent ways to achieve the desired outcome with highest 
probability, measured by average (equilibrium) or MFPT 
(noncquilibrium) attributes. 
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